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Introduction - The oscillator basis is widely used in nu- 
clear structure computations because it allows the prac- 
titioner to exploit and implement all symmetries of the 
nuclear many-body problem, and because its localized 
nature corresponds well to the structure of the self-bound 
atomic nucleus. After all, the nuclear shell model is based 
on the harmonic oscillator with a strong spin-orbit split- 
ting [T]. Several computational implementations of ab 
initio methods [5J [3] and nuclear density functional the- 
ory [4j essentially start from the oscillator basis. Ba- 
sic observables sought in such computations include the 
binding energies and radii. Ideally, the computed ob- 
servables should be independent of the parameters of the 
employed oscillator space, i.e., the maximum number of 
oscillator quanta N and the frequency of the oscillator 
wave functions. This ideal is often difficult to reach in 
practice, so various empirical extrapolation schemes [S]-[5] 
have been applied, but all lack a firm theoretical founda- 
tion. 

The proper accounting for corrections to nuclear ener- 
gies and radii that arise in finite oscillator spaces is an 
important problem for several reasons. First, a theoreti- 
cal foundation of such corrections would enable the prac- 
titioner to extrapolate reliably from smaller model spaces 
and thus extend the reach of some computational meth- 
ods. This is particularly important for weakly bound 
nuclei where the Gaussian falloff of the oscillator ba- 
sis can capture a halo state often only in unachievable 
large model spaces [B]. Second, uncertainty quantifica- 
tion of results - standard in experimental research - is 
increasingly taking place in nuclear structure theory |10j . 
Here, the quantification of theoretical uncertainties due 
to the nuclear interaction is possible for interactions from 
effective field theory (EFT), but the robust quantifica- 
tion of errors due to finite oscillator spaces is lacking. 
Finally, important steps towards an harmonic-oscillator 
based EFT for the nuclear shell model have been made 
recently [TTI - fT5] . Such a theory should also control and 
exploit the limitations of the finite model space. 

In this Rapid Communication, we derive corrections of 
nuclear energies and radii that are due to finite oscillator 
spaces. We build on the insights of Coon et al. [TT], who 
focus on the infrared and ultraviolet cutoffs induced by 
a truncated basis. Our derivations are based on simple 
arguments and verified in a one-dimensional model. We 
apply the results to -'^^O and demonstrate that the the- 



oretical corrections agree well with the numerical data. 
Calculations for the ^He ground-state energy and neu- 
tron radius show that predictions are feasible even for 
halo nuclei. 

Theoretical derivation - For a particle in a box with 
periodic boundary conditions, Liischer derived the cor- 
rections to bound states due to the finite size of the 
box [13]. Our derivation is analogous, except that the 
size of the box is now given in terms of the spatial exten- 
sion of the oscillator basis and we deal essentially with 
Dirichlet boundary conditions. Let us consider a model 
space of oscillator wave functions with maximum oscil- 
lator energy E = hn{N + 3/2). In practice, one has to 
choose hil and N such that the momentum cutoff A of 
the employed interaction is smaller than the ultraviolet 
(UV) momentum 

Auv = V^iN + 3/2)h/b , (1) 

and that the radius r of the nucleus is smaller than the 
radial extent 

Lo = V2(A^ + 3/2)6 (2) 

of the employed oscillator space. Here, b = ^ h/{mfl) 
is the oscillator length of our basis, and m denotes the 
nucleon mass. The definitions ([I]) and ([2| are indeed the 
maximum momentum and displacement, respectively, of 
a particle in a harmonic oscillator at energy E ~ hH{N + 
3/2). They differ from previous scaling relations [T51 116] 
by factors of \/2. 

In practice, satisfying the UV condition A < Auv and 
the infrared (IR) condition r < Lq does not guarantee 
converged nuclear structure results in the oscillator ba- 
sis because the momentum cutoff A is usually not sharp, 
and the nuclear wave function extends beyond the nuclear 
radius r. However, nuclear interactions from chiral EFT 
and from renormalization group transformations exhibit 
a super-Gaussian falloff in momentum space, whereas the 
nuclear wave function only falls off exponentially in co- 
ordinate space. Thus, once A < Ajjv holds, the UV 
convergence in momentum space will be rapid, and one 
is dominated by corrections from the slower falloff in co- 
ordinate space. Practitioners of nuclear structure com- 
putations know this very well (see, e.g.. Fig. |6] below): 
When energies are plotted as a function of hfl, the min- 
imum initially shifts toward larger values of hfl as N 



2 



is increased. However, once UV convergence has been 
reached, further increasing TV shifts the minimum back 
to lower values of HQ to capture the coordinate-space tail 
of the wave function. In what follows, we will assume UV 
convergence and compute the correction from incomplete 
IR convergence. 

The finite extent of the oscillator basis up to a radius L 
in coordinate space essentially requires the wave function 
to vanish at r w L. The maximum radius Lq from Eq. ([2| 
is only a leading-order (or asymptotically valid) estimate 
because the oscillator wave function decays rapidly be- 
yond the classical turning point. An improved estimate 
for L using the intercept of the tangent at r = Lq is [13 



L w Lo +0.54437 &(Lo/&) 



-1/3 



(3) 



which we use with the analogous expression for Ajjv in 
numerical examples below. For our derivation of IR cor- 
rections, we adapt the discussion in Ref. [18j . Given a 
boundary condition at r = L beyond the range of the 
nuclear potential, we write the energy compared to that 
for L = oo as 



El — En 



AEr 



(4) 



and we seek an estimate for AE^, which is assumed to 
be small. 

Let UE{r) be the radial solution with regular boundary 
condition at the origin and energy E. We denote the par- 
ticular solutions UE,^{r) = UL{r) and UE^{r) = Uaoif). 
Then the linear energy approximation is (for r < L) |18| 



UL{r) w Uoo('') + AEl 



duE{r) 



dE 



(5) 



assuming a smooth expansion oi ue about E — E^o at 
fixed r. Evaluating Eq. ([s]) at r = L with the boundary 
condition ui^{L) — 0, we find 



Substituting Eq. ([o]) at r = L into Eq. we obtain 



AEj 



dOLE 




dE 





The prefactor in the square brackets depends on details of 
the interaction (but not on L), and will be fit to numerical 
data when Eq. (10) is used together with Eq. Q. Thus, 
the main result is 



El — En 



floe 



(11) 



and in practical applications one can treat Eoo^ clq and 
fcoo (in cases where the separation energy is not known) 
as fit parameters. Note that our result (11) explains the 
exponential decay observed empirically in Ref. |llj . In 
contrast to the Liischer result in which the energy is al- 
ways lowered by periodic images of the potential [14j . 
the energy from Eq. ( 11 ) is always increased by the shift 
of a node from r = oo to r = L, consistent with the 
variational nature of the truncated basis expansion. 

Let us next turn to radii. It is convenient to express 
the radius squared as the infinitc-model-space result plus 
a correction term 



{r')L 



A(r^) 



where 



2n ^ Io \uL{r)\^ dr \uno{r)\^ r"^ dr 



(12) 



(13) 



Because the dependence on L of UL{r) in Eq.(5) is con- 
fined to AEl, when is substituted into Eq. (13) the 



L dependence of each separate integrand comes entirely 
from the upper integration limit. Therefore we can use 
the asymptotic expressions Uoo(r) — > AooC"*^""'' and 



AEr 



-Uoo[L) 



duE{L) 



dE 



(6) 



duE{r) 
dE 



A. 



AEr 



(14) 



which is the estimate we seek. For general E, the asymp- 
totic form of the radial wave function for r greater than 
the range R of the potential is 



UE\r) — > AE{e 



- aEG 



+kEr 



(7) 



with the known case Uoo (r) A^e-''^'- for E = E^o- 
Here, fcoo is determined by the nucleon separation energy 



^ 2m 



(8) 



We take the derivative of Eq. ^ with respect to energy, 
evaluate at E — Erx, using aE^ — and dkE/dE — 
—m/{'h?kE) to find 



duE{r) 



dE 



daf 



dE 



0\ 



(9) 



to identify the leading-order expression A(r^)i (x 
(r^)oo(2fcooi)^e~^''°°^. (Note that any L-independent 
terms are guaranteed to cancel by the definition 



Eq. (13).) The next-to-leading-order expression scales as 
(2fcooi) exp (— 2fcoo-^) because the condition ul{L) = 
ensures there is no quadratic term in 2kooL. Thus, the L 
dependence of the squared radius is (with /3 = 2kooL) [T7] 



Here, 



Co, 



r')oo[l-(co/?' + Ci/3)e-^] . (15) 
and ci are fit parameters while koo 



should be determined in fitting the energy ( 11 ). The ap- 



proximation ( 15 ) is valid in the asymptotic regime /3 3> 1. 

> 3 because the leading-order 
3, and the next-to- 



In practice, one needs 
correction has its maximum at /3 
leading order corrections is approximately suppressed by 
one order of magnitude for /? ^ 3 (with cq and ci of order 
unity) . 



3 



Equations (11) and (15) are the main results of this 
Rapid Communication. A few comments are in order 
before we turn to appUcations of these results. Note that 
we derived these results in the laboratory system. For the 
nuclear ^-body problem, we could also have exploited 
the separation of the center-of-mass coordinate in the 
oscillator basis and followed a similar derivation for the 
A*'^ particle with respect to the center of mass of the 
remaining (A—l) particles. This would rescale L and the 
momentum of the A^^ particle accordingly, but the final 
results are unchanged when re-expressed in laboratory 
coordinates. In light of the approximations involved in 
defining L and the corrections to the energy and radius, 
in actual fits to numerical results one might want to treat 
/coo StS a fit parameter even when the separation energy 
is known. 

Applications - As a first check, we consider a toy 
model in one dimension with the Hamiltonian H = 
— wocxp (— x^). Here, x is given in units of the os- 
cillator length b. In one dimension, the constant 3/2 in 
Eqs. ([T]) and ([2| has to be replaced by 1/2. We com- 
pute the ground-state energy and the squared radius for 
?;o = 0.5 in large oscillator spaces to obtain fully con- 
verged results for the ground-state energy and the radius. 
In this simple case, the ground-state energy is given in 
terms of the separation energy ([s]) as i?oo = —S. Fig- 
ure l] shows the correction A(r^) as a function of L. The 
dashed line results from a leading-order fit to Eq. (15), 



and the agreement between numerical data and the the- 
oretical prediction extends over ten orders of magnitude. 
The inset shows that the L-dependcnt energy correction 
also agrees with the prediction (11). 




L/b 



FIG. 1: (color online) Finite-basis-size correction of the 
squared radius (crosses) compared to Eq. ( |15[) (dashed line) 
for a toy model. The squared radius is r^/6^ ~ 1.736, and 
koob ~ 0.595. Inset: Finite-basi s-siz e correction of the energy 
(data points) compared to Eq. (10) (dashed line). 



Let us turn to the nuclear many-body problem. We 
employ the nucleon-nucleon interaction from chiral EFT 



by Entem and Machleidt (TH], and compute the ground- 
state energy and radius of the nucleus "'^^O with the 
coupled-cluster method in its singles and doubles ap- 
proximation with triples corrections |15l I20j . We employ 
model spaces with frequencies 42 < h^l/MeV < 76 and 
with N = 12, 14. To ensure that the computed results 
are practically UV converged, we only use those oscilla- 
tor spaces for which Auv is sufficiently large. Figure [2] 
shows the results for the ground-state energy as a func- 
tion of L. The circles, up triangles, and down triangles 
denote points with Auv > HOO MeV, Auv > 1200 MeV, 
and Auv > 1300 MeV, respectively. The points all fall 
on a line because UV convergence has practically been 
achieved. Thus, we can apply our theory. The lines show 
fits to Eq. ( 11 ) with fit parameters -Eqo: s-^d fcoo- Note 
that the result Eao ~ —122.6 MeV of the fit depends very 
weakly on Auv, the difference being about 0.2 MeV. In 
the fits, we obtain fcoo ~ 0.95 fm~^, and this agrees well 
with the decay of the Pi/2 orbital that contributes to the 
density 21j. 
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FIG. 2: (color onhne) Open symbols: Ground-state energy of 
as a function of L. Lines: Fits to Eq. Ill) yield Eoo ~ 
fm-i. ■ " 



-122.6 MeV and fc„ 



0.95 



Auv from Eq. fll 



Next we consider the radius. We use Eq. ( 15 ) including 



the next-to-leading order correction and fit the parame- 
ters {r'^)oo, Co, and ci to data, with fcoo taken from the 
fit of the ground-state energy. The result is shown in 
Fig. [3j The circles, up triangles, and down triangles de- 
note points with Auv > HOO MeV, Auv > 1200 MeV, 
and Auv > 1300 MeV, respectively. The lines show the 
corresponding fits and asymptotes, and the extrapolated 
radius is r sa 2.34 fm. It is particularly satisfying that 
the extrapolation also works well for the few data points 
with Auv > 1300 MeV. 

We also consider the challenging case of a halo nucleus. 
The isotope ^He is only bound by about 0.97 MeV with 
respect to ^He, and thus exhibits a two-neutron halo. 
Note that ^He is not a bound nucleus, and that the neu- 
tron separation energy of ^He is about 1.86 MeV. As a 
consequence of the weak binding, the matter radius of 
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FIG. 3: (color online) Open symbols: Squared radius as a 
function of L for ^^O. Lines: Fits of Eq. ( |l5[ ) with kac fixed 
from the energy fit. A.uv from Eq. ([T]). 



^He is unusually large (about 2.4 fni compared to 1.5 fm 
for *He) [15H21|. We address this challenge by applying 



the finite-basis-size corrections to the energy (11) and 
neutron radius (15). 



Our test case uses NCFC results [B] obtained for a 
chiral EFT nucleon-nucleon interaction that was soft- 
ened via a similarity renormalization group (SRG) trans- 
formation with a parameter A — 2.0 fm""'^. Fig- 
ure |4] shows the fit of the ground-state energy for model 
spaces with Kuv > 660 MeV and fifi > 24 MeV 
(which ensures a small UV correction). The fit yields 
~ —29.87 MeV, and the computed two-neutron 
separation energy is about 0.95 MeV. Thus, both en- 
ergies are in good agreement with experiment (despite 
the absence of a three-body force). The fit also yields 
hk^ « 93 MeV corresponding to a neutron-separation 
energy of S* sa 4.6 MeV, which is significant larger than 
the experimental value. Of course, the unbound nucleus 
^He cannot be computed reliably in the oscillator basis, 
but the interpretation of in this case requires further 
study. 

Now we can extrapolate the neutron radius of ^He. 
The results of NCFC calculations are shown in Fig. [5] 
Without a knowledge of the finite-basis-size corrections, 
it would be impossible to make any reasonable prediction 
for the radius because of the apparent lack of conver- 
gence. Note, however, that we are in the UV-converged 
regime and 2kooL > 3 for the NCFC data points in Fig.|4] 
Thus, we should be able to apply our correction formula. 
The solid and dashed lines in Fig. [5] show the results for 
the radius based on a fit at leading order and next-to- 
leading order, respectively. At next-to-leading order we 
find r « 2.40 fm, and this prediction is in reasonable 
agreement with deductions from data |23j . Several ad- 
ditional points at large L not included in the fit are in 
good agreement with the extrapolation. 

Our theoretical results have been derived under the 
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FIG. 4: (color online) Circles: NCFC ground-state energies of 
the halo nucleus ''He from Ref. [6j. Full line: Fit of Eq. jll]) 
yields ^ -29.87 MeV (dashed line) and hkac ^ 93 MeV. 
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FIG. 5: (color online) Circles a nd squares: NCFC neutron 
radii of ®He. Lines: Fits of Eq. (151 to circles with fcoo fixed 



from energy fit yield a radius of 2.37 fm at LO and 2.40 fm at 
NLO. 



assumption that UV convergence has been reached. It 
would also be useful to know finite-basis-size corrections 
in the opposite regime where IR convergence is estab- 
lished, for instance through calculations in model spaces 
with sufficiently small hfl. The remaining UV correc- 
tions would, of course, depend on the interaction at hand. 
We have not yet established a theoretical derivation but 
can resort to empirical findings for SRG-evolved nucleon- 
nucleon interactions from chiral EFT with an SRG pa- 
rameter A j6j. We plotted IR converged ground-state 
energies (computed at large values of 2kaoL) for vari- 
ous light nuclei as a function of Kjjv and found as in 
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Refs. [11] [12] that the empirical formula 



E{Auv) = Eoo+Aoe' 



(16) 



works quite well [T^. (In practice we allow A to be a fit 
parameter to optimize the fit.) This formula is consistent 
with an empirically successful ansatz used for individual 
hn values (e.g., see Ref. [B]). However, Eq. (161 allows 
results with different hil. to be fit all at once. 
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FIG. 6: (color online) Data points connected by full lines: 
NCFC ground-state energies of ^He. Red dashed line: Fit 
of Eq. ([17| to NCFC data with iV^ax = 6-10 yields « 
—29.84 MeV. Data points connected with dotted line: Ap- 
plying the correction of Eq. ([T7| to the NCFC data. (iV = 

iVmax + 1) 



If we combine the empirical UV formula and the the- 
oretically founded IR formula assuming the corrections 
are approximately independent, then 

E{Auv, L)^E^+ Aoe-^A^'v-M? + A^e-^"-'' . (17) 

Here i?oo, ^ij and fcoo are fit parameters that are 
determined from a simultaneous optimization to data at 



all Ml, including in the intermediate region where both 
IR and UV corrections are significant. The resulting 
value « -29. 84 MeV from using aU A^,„ax = 6-10 
points, which is in good agreement with the IR-only fit 
in Fig. [4| is plotted as a dashed red line in Fig. [6| The 
points connected by dashed lines are obtained by sub- 
tracting the corrections in Eq. (17 1 from the NCFC ener- 



gies. Thus a perfect fit would find all points lying on the 
line for Eoo- (Note: the A^max values in the figure are for 
excitations above the ground state [6] , so TV = iVmax + 1 
for ^He [IT1[T7|.) All corrected points included in the fit 
are close to the Eoo line and even the corrected Nj^ax = 4 
energies (which were not included in the fit) are only 
slightly overbound. 

In summary, we derived analytical results for the finite- 
basis-size corrections of nuclear radii and energies that 
are valid in oscillator spaces with converged ultraviolet 
physics. The computation of the corrections is robust 
and appears to be applicable to halo nuclei. In combi- 
nation with an empirical formula for the ultraviolet cor- 
rection for SRG-transformed interactions, consistent and 
much-improved extrapolations of ground-state energies 
are possible. The analytical results can be extended to 
other long-range observables that are sensitive to the tail 
of the nuclear density. A systematic study of the ex- 
trapolation procedure including an error analysis will be 
presented in a future work |17| . 
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